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Based on a true phase space probability distribution function and an ensemble averaging pro- 
j^jf)' cedure we have recently developed [ Phys. Rev. E 65, 021109 (2002) ] a non-Markovian quantum 

r-{ i Kramers' equation to derive the quantum rate coefficient for barrier crossing due to thermal ac- 

■ tivation and tunneling in the intermediate to strong friction regime. We complement and extend 

this approach to weak friction regime to derive quantum Kramers' equation in energy space and 
the rate of decay from a metastable well. The theory is valid for arbitrary temperature and noise 
correlation. We show that depending on the nature of the potential there may be a net reduction 
(— I | of the total quantum rate below its corresponding classical value which is in conformity with earlier 

. observation. The method is independent of path integral approaches and takes care of quantum 

D ' effects to all orders. 
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I. INTRODUCTION 



The dynamics of noise-induced rate processes was first successfully treated in a seminal paper by Kramers in 
\ 1940 0]. With the advances in experimental methods for monitoring ultrafast processes on microscopic spatial and 
temporal scales over the last two decades this has been the subject of numerous investigation from classical, 

semiclassical and quantum mechanical point of view B-R]. The classical Kramers' th eory has thus been extended 
to non-Markovian dissipation models [pH10||, generalization to complex potential [pd|- [l3| and to many degrees of 
freedom |L4|,[l5), fluctuating barrier problem and non-stationary activated processes ]17|], thermal ratchet |L8 19 
and molecular motors [f20f, a nalysis of quantum~P.|5| ,pT|-p3| and semiclassical effects |24| , calculation of time-dependent 
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transmission coefficient J25| , P6|1 , fractional kinetics (g7fl28[^nonequilibriurn open systems p9|.p0|] , activationless escape 
of a free Brownian particle jpl) and other related issues 

Although classical Kramers' equation was proposed more than sixty years ago and quantum Kramers' problem 
of escape from a metastable state has attracted wide attention over the last two decades Q, the quantum version 
of Kramers' equation was not reported in the literature. This is probably because of the fact that the traditional 
^vq ' method of treatment of quantum Kramers' problem rests on calculation of partition function for a system-reservoir 
Hamiltonian in terms of path integrals, rather than on evolution of probability distribution function as used in classical 
theory of stochastic processes. Very recently we have developed ]3^ | a method based on true quantum (c-number) 
phase space distribution function (rather than quasi-probability function, like Wigner function |33| ) to derive for the 
first time an exact non-Markovian quantum Kramers' equation which is valid for arbitrary temperature and friction. 
The solution of this equation as an appropriate boundary value problem results in an expression for quantum rate 
coefficient which not only reduces to Kramers '-Grote-Hynes rate in the classical limit but also to the result 
corresponding to zero-temperature tunneling in the full quantum limit, treated by Caldeira and Leggett |34|] in early 
, eighties. The rate coefficient thus derived pertains to spatially-diffusion limited processes and is therefore valid for 
intermediate to strong friction regime. We undertake the present study with the following specific objectives to 
. ,— i | complement this work in the low friction regime where the process is controlled by energy diffusion. 

(1) To extend the treatment of quantum Kramers' problem for low to low-moderate friction we develop a quantum 
Kramers' equation for energy diffusion which is a quantum version of classical non-Markovian equation of Carmeli 
and Nitzan (CN) |n]] proposed in early eighties. 

(2) Our aim here is the inclusion of memory effects for arbitrary noise correlation of the heat bath kept at an arbitrary 
temperature taking into consideration of the quantum effects (corrections) to all orders. 
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(3) We solve the quantum Kramers' equation for energy diffusion to derive an explicit form of rate coefficient in the 
weak friction regime and show that it reduces to non-Markovian counterpart of Hanggi and Weiss |35[] in the classical 
limit. Furthermore it provides the rate coefficient of low temperature tunneling (down to absolute zero) in the quantum 
limit. The present theory thus interpolates between thermal activation and tunneling for weak dissipation within a 
single scheme and is a direct extension of classical theory to quantum domain. 

The classical non-Markovian Fokker-Planck equation in the energy variable for arbitrary noise correlation was first 
proposed by CN fllPf . The detailed classical analysis by several groups ]Tc| , |35| -|37t revealed that the rate, in general, 
is significantly modified by memory effects when compared to corresponding Kramers' theory in the static friction 
limit. As mentioned earlier, the traditional quantum treatment of the Kramers' problem in weak friction limit on 
the other hand, is based on functional integral approach ||] which takes care of dissipative tunneling |34|]. Since 
for weak friction limit at a finite temperature one finds a small population at the upper energy levels of the system 
which results in non-equilibrium effects, quantum correction to classical Kramers' weak damping results above the 
cross-over temperature is of considerable interest. Several authors p8|-|43]| have addressed this problem in relation 
to nonequilibrium quantum tunneling out of a metastable state. Although the method of functional integrals as 
employed in these theories has been successful in treating arbitrary coupling and correlation time scales in a formally 
exact manner, analytic evaluation of these integrals usually requires semi-classical approximations, e.g., semiclassical 
steepest descent method, WKB approach, etc., or other specific cases, which put restriction on applicability of the 
theories in several ways. Secondly, the weak coupling theories ( which have been extensively used in quantum optics 
jfJJ since seventies ) based on master equations that make use of quasi-probability functions, like Wigner function |3^| 
often pose serious difficulties concerning negativity or singularity of the probability distribution functions as discussed 
in detail in earlier work p5j |. Third, when the system or the system-reservoir coupling is nonlinear, the differential 
equations concern higher (than second) derivatives of quasi-probability functions p6| ] for which quantum-classical 
correspondence gets blurred. Our approach here is based on true quantum probability phase space function and is free 
from such difficulties. Furthermore an important decisive advantage of the scheme is that it allows us to implement 
the classical non-Markov theories of activated processes in a full quantum setting without taking recourse to any 
semiclassical technique. In what follows we specifically apply the classical procedure of Lax 47 1, CN jio]], Biittiker, 
Harris and Landauer (BHL) fl48|], Hanggi and Weiss p5| in our quantum phase space formulation to develop a non- 
Markovian quantum Kramers' equation in energy variable and derive an expression for quantum rate coefficient in 
the spirit of classical Kramers' theory. The quantum Kramers' equation and the rate coefficient are classical looking 
in form but quantum mechanical in their content and it is easy to recover their non-Markovian classical counterparts 
in the limit fi — > 0. 

The Kramers' kinetics in the low friction regime is just not a theoretical issue today but has been a subject of 
experimental investigation over the last two decades p|,[49|-|53"[ . A number of experimental work in chemistry aimed 
at detecting Kramers' turnover phenomena, in various reactions which can be conveniently explained in terms of a 
one-dimensional model, e.g., iodine atom recombination in various inert solvents p9[ , chair-chair isomerization of 
cyclohexane |30| , excited state isomerization of 2-alkylanthracene |5^] . Another class of experiments where the energy 
diffusion mechanism has been successfully implemented concern resonant activation of a Josephson junction p7|fl and 
decay of zero voltage state in a current-biased Josephson junction ]5^| . The non-exponential decay behaviour in spin 
glass [U is also an area of active research in this context. All these problems have their quantum counterparts which 
are being considered for further studies in rate theory although the experimental evidence of some of the theoretical 
predictions in low temperature quantum effects in weak friction regime is still awaited. 

The outlay of the paper is as follows. We introduce a c-number representation of a g eneralized quantum Langevin 
equation in Sec. II. This formulation helps us to use the classical formulation of CN |10| for deriving a non-Markovian 
quantum Kramers' equation in energy space in Sec. III. We solve the problem of quantum energy diffusion controlled 
rate coefficient in the spirit of classical theory |}5 48 in Sec. IV. This reduces to classical rate expression of Kramers- 
Hanggi- Weiss p5| form in the limit Ti — > 0. An explicit example with a cubic potential is worked out to illustrate the 
theory in Sec. V. The paper is concluded in Sec. VI. 



II. THE QUANTUM GENERALIZED LANGEVIN EQUATION (QGLE) IN C-NUMBERS 

We consider a particle in a medium. The latter is modeled as a set of harmonic oscillators with frequency {coi}. 
Evolution of such a quantum open system has been studied over the last several decades under a variety of reasonable 
assumptions. Specifically our interest here is to develop an exact description of quantum Brownian motion within the 
perview of this model described by the following Hamiltonian p4| , 
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H=— + V(X)+J2 
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(2.1) 



Here X and P are co-ordinate and momentum operators of the particle and the set {qj, pj} is the set of co-ordinate 
and momentum operators for the reservoir oscillators coupled linearly to the system through their coupling coefficients 
Kj. The potential V(X) is due to the external force field for the Brownian particle. The co-ordinate and momentum 
operators follow the usual commutation relation [X,P] = ifr and [qj,pj] = ihdij. Note that in writing down the 
Hamiltonian no rotating wave approximation has been used. 

Eliminating the reservoir degrees of freedom in the usual way p4|j5^j5^ ] we obtain the operator Langevin equation 
for the particle, 



X(t)+ / dt'(3{t-t')X{t') + V'{X)=F{t) 



(2.2) 



where the noise operator F(t) and the memory kernel /3(t) are given by 



F(t) = [{?j(0) - X(0)\ kj cosujjt + pjiO)*, 1 / 2 shiLUjt 



(2.3) 



and 
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(2.4) 



with Kj = Wj ( masses have been assumed to be unity ). 

The Eq.(2.2) is an exact quantized operator Langevin equation which is now a standard textbook material 0] 
and for which the noise properties of F(t) can be defined using a suitable initial canonical distribution of the bath 
co-ordinates and momenta. Our aim here is to replace it by an equivalent QGLE in c-numbcrs. Again this is not a 
new problem so long as one is restricted to standard quasi-probabilistic methods using, for example, Wigner functions 
|j33f . To address the problem of quantum non-Markovian dynamics in terms of a true probabilistic description we, 
however, follow a different procedure. We first carry out the quantum mechanical average of Eq.(2.2) 
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dt'(3(t-t')(X(t')) + (V'(X)) = (F(t)) 



(2.5) 



where the average (. . .) is taken over the initial product separable quantum states of the particle and the bath oscillators 
at t = 0, |</>){|ai)|a2) ■ ■ • |ckjv}}- Here \<j>) denotes any arbitrary initial state of the particle and |a,) corresponds to 
the initial coherent state of the i-th bath oscillator. |a,) is given by |a,-) = exp(— |ai| 2 /2) ^^ =0 (a™ i /v / "^)l n i)j 
oti being expressed in terms of the mean values of the co-ordinate and momentum of the i-th oscillator, (<ji(0)) = 
(Vh/2coi)(ai + a*) and (pi(0)) = i^/hujjj2(a* — oti), respectively. It is important to note that (F(t)) of Eq.(^) 
is a classical-like noise term which, in general, is a non-zero number because of the quantum mechanical averaging 
over the co-ordinate and momentum operators of the bath oscillators with respect to the initial coherent states and 
arbitrary initial state of the particle and is given by 
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(2.6) 



It is convenient to rewrite the c- number equation ( 2^5 ) as follows 



(X(t)) + / df(3(t - t')(X{t')) + (V'{X)) = F(t) 
Jo 



(2.7) 



where we let the quantum mechanical mean value (F(t)) — F(t). We now turn to the second average. To realize F(t) 
as an effective c- number noise we now assume that the momenta (pj(0)) and the shifted co-ordinates {(<?j(0)) — (X(0))} 
of the bath oscillators are distributed according to a canonical distribution of Gaussian forms as 



Vj = AAexp 



2huij (rij 



(2.8) 
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so that for any quantum mechanical mean value Oj((pj(0)) , {(<?j(0)) — (X(Q))}) the statistical average (. . .)s is 



rij indicates 



(Oj)s = J O,(fe(0)),{(%(0)>-(1(0))}) 

x^(fe(0)>,{fe(0))-(l(0))}) 
xd<fc(0))d{(&(0))-(X(0))} . 

average thermal photon number of the j-th oscillator at temperature T and 



Here 

l/[exp {hujj /ksT) — 11 and M is the normalization constant. 

The distribution (UM) and the definition of statistical average (2.9) imply that F(t) must satisfy 



(F(t))s = 



and 



(F(t)F(t')) s = ^KjJkjj fcoth^) 



COS U)j(t — t') 



(2.9) 



(2.10) 



(2.11) 



That is, the c-number noise F(t) is such that it is zero centered and satisfies the standard quantum fluctuation- 
dissipation relatio n (F DR) as known in the literature |55) in terms of quantum statistical average of the noise operators. 
The distribution ( ]2.8| ) is thus an ansatz introduced to calculate the ensemble a vera ge over the quantum-mechanical 
mean values of the bath oscillators. Its justification lies in the fact that with ( |2.9| ) it reproduces the correct noise 
properties of the bath, i.e., the quantum fluctuation-dissipation relation for the c-number quantum noise ( 2.1 1| ) along 
with ( [2.10| ). Secondly, the distribution (2.£) has a form which is Boltzmann-like (but not a Boltzmann distribution) 
since the width parameter of the Boltzmann distribution kT gets replaced by 7kj ,(ri,- + §) 

To proceed further we now add the force term V'((x)) on both sides of Eq.(2.7) and rearrange it to obtain formally 



x(t) + / dt'/3(t - t')x(t') + V'(x) = F(t) + Q{x, t) 



(2.12) 



where we let (X(t)) = x(t) for simple notational convenience and 

Q(x,t) = V'(x)-(V'(X)) 



(2.13) 



represents the quantum mechanical dispersion of the force operator V'(X) due to the system degree of freedom. 
Since Q(x,t) is a quantum fluctuation term Eq.( 2.12 ) offers a simple interpretation. This implies that the classical 
looking QGLE is governed by a c-nu mber quant um n oise F(t) which originates from the quantum mechanical heat 



bath characterized by the properties ( 2.1C ) and ( 2.11 ) and a quantum fluctuation term Q(x,t) due to the qua ntum 
nature of the system characteristic of the nonlinearity of the potential. Although because of the last term in Eq.(2.12) 
the equation looks formal and implicit, the actual structure of Q(x,t) gets more transparent as we go over to the 
beginning of the next section. A recipe for calculation of Q(x,t) is given in Ref. 45 57|j5^| 



We summarize the above discussions to point out that it is possible to formulate an exact QGLE ( 2.12 ) of the 
quantum mech anical mea n value of position of a particle in a medium, provided the classical-like noise term F(t) 
satisfies (2.10) and ( 2.1l| ). The important new content of the approach is that to realize F(t) as a noise term 
we have split up the standard quantum statistical averaging procedure into a quantum mechanical mean (...) by 
explicitly using an initial coherent state representation of the bath oscillators and then a statistical average (. . .)s of 
the quantum mechanical mean values with distribution ( |2 .8[ ) . This is distinctly different from the usual procedure 
of quantum statistical averaging where the quantum mechanical average is carried out with number states over the 
noise operators followed by an ensemble average with Boltzmann distribution. Two pertinent points are to be n oted : 
First, it may be easily verified that the distribution of quantum mechanical mean values of the bath oscillators (2.5) 



reduces to classic al M axwell-Boltzmann distribution in the thermal limit, Tiujj <C ksT. Second, the vacuum term in 
the distribution (2.8) prevents the distribution of quantum mechanical mean values from being singular at T = 0; 
or in other words the width of distribution remains finite even at absolute zero, which is a simple consequence of 
uncertainty principle. The procedure has been recently implemented by us to formulate a quantum theory of Brownian 
motion [ fig ] and to propose an exact non-Markovian quantum Kramers' equation |32j with true probability distribution 
functions. 
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III. QUANTUM KRAMERS' EQUATION IN ENERGY SPACE 



Let us begin by noting that the generalized quantum Langevin equation ( 2.12| ) of a Brownian particle in presence 
of an external force field takes into account of arbitrary coupling between the system and heat bath and contains 
quantum corrections, Q(x,t) due to system to all orders. To make the later assertion explicit we now express the 
operators X and P as 



X(t) = (X(t))+SX(t) 
P(t) = (P(t))+6P(t) 



(3.1a) 



By construction (SX(t)} = 0, (SP(t)) = and [SX,8P] = ih. Expanding (V'(X)) around (X) (= x) in a Taylor series 
we obtain 



1. 



(V'(X)) = V'{x) + -V"'(x)(6X 2 ) + 



(3.1b) 



Therefore Q(t) can be expressed as 



oc 1 

= "E j^jy V n (x)(5X( n -V(t)) 



(3.1c) 



Here V n (x) denotes the n-th derivative of the classical potential. The role of Q(x,t) is therefore to modify the classical 



potential V{x) in Eq.(2.12). Q(x, t) can be calculated by solving (5X n (t)) order by order. To the lowest order (second) 
(X) and (SX 2 ) follow a coupled set of equations as given in Eqs.55(a-e) of Ref. f45|] . (Higher order equations, e.g ., the 
fourth order equations, are given in Ref. 58 1). For convenience, we will now split up the right hand side of Eq.( 3.1c ) 
into a time-independent and a time-dependent part as, 



oo 

0(*) = -E? W y n (x)(SX^(0)) + g(t) 



(3. Id) 



where 



n=3 >■ >' 



For future use it is convenient to write g(t) in the Taylor series of the form, 



(3.1e) 



1 t" 



OO OO 

= " EE (n _i)!, ;i! 

m=0n=3 v ' 



-V 



(3.1f) 



The Langevin equation (2.12) then reduces to 



Expressing 



+ f drP(t - r)i(r) + V'{x) + £ —!—V n (x){dX^-V(0)) = F(t) + g(t) 



oo 1 

V q (x) = V(x) + J2 j—fyV^WiSX^iO)) 



Eq.(^2[) takes the form 



n— 3 



(3.2) 



(3.3) 



(3.4) 



v + / dr0(i - t)v(t) + Vfa) = F(t) + g{t) 
Jo 



(3.5a) 
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Eg. (||a|) is our starting Langevin equation. The potential V q (x) appearing in ( |3.3| ) and ( |3.5a| ) is not the classical 
potential but a renormalized on e wit h quantum corrections. The damping memory kernel (p~4 ) is identified by the 
fluctuation-dissipation relation ( [2.11 ) by noting that in the continuum limit 



(F(t)F(0)) s = 



dojK(u>)p(u>)hu! ( coth 



and 



p{t) 



huj 



duj k{lo) p{(jj) cosujt 



cosuut = C(t) 



in the Fourier domain can be related as 



C c (uj) = \hu (coth 



TlLU 

2k B T 



(3.5b) 



(3.5c) 



(3.5d) 



where C c (uj) and /3 c (w) are the cosine transforms of C(t) and /3(i), respectively. For convenience we now express the 
Fourier transform of (3(t) as 



dt(3{t) exp(— imut) 



(3.6) 



We now consider the following time scales in the dynamics relevant for energy diffusion in the weak friction limit, 

7 < 1/t c < u (3.7) 

where 7 is the friction arising due to interaction with the bath, evaluated in the Markovian limit. r c is the correlation 
time of the noise due to heat bath and to is the li near ized system frequency, which for a Brownian particle is assumed 
t o be very high. This separation of time scales in (3/7) and casting of an operator Langevin equation in c- number form 
(3.4-3. 5a) allow us to implement a classical met hod for s olving the problem of quantum energy diffusion. Following 
the standard procedure one can transform Eqs.( 3~^ - 3.5a ) to the action (J) and angle (0) co-ordinates with the help 
of a Jacobian matrix as 



J 



_ dv dx 

d<p d$ 

dv dx 

dJ dJ 

dv dx 



O.J 



dJ 



OH 



ox-fidT/3(t-T)v(T)+F(t) + g(t) 



Thus we have 



J 



O.r 



(j) = uj{J) - 



dr(3(t - t)v{t) + F(t) + g(t) 
dx r '* 



dJ 



drf3{t - t)v{t) + F(t) + g{t) 



(3.8) 
(3.9) 



Here v represents the velocity of the particle. For the deterministic part of the system's Hamiltonian given by 
H = [l/2)v 2 + V q (x) we may write, 



dH{J) 
dJ 



(3.10) 



Sinc e V q {x) (see Eq.( p.3[ )) contains quantum corrections, our J and <fi are quantum (c- number) variables as implied in 
(3.3). In the absence of quantum corrections they become classical variables of CN p0| . The canonical transformation 
from (x,v) space to (J, <fi) space has been done with the deterministic Hamiltonian. We can therefore expand x and 
v as, 



x(J, cf>) 
v(J,<j>) 



2J x n (J) exp(in^) 

n— — 00 
00 

^2 v n (J) exp(incf)) 



(3.11a) 
(3.11b) 



n— — 00 



G 



with 



x n = x* and v n = v* 



(3.12) 



Differentiating Eq.(3.11a) with respect to time and noting that in the action-angle variable space cj> = lu(J) we can 
write, 



,(J) = inu(J)x n (J) 



(3.13) 



Since, we are considering the motion in one dimension only, we can choose J and 4> m such a way that we can make 
the simplification for x — x as, 



- [x n exp(m^) + x* L exp(— inq 



n— — oo 



Inserting Eq.(3.12) we get, 



^ :v 

= — [x n exp(in<j>) + x- n exp(—inq 



n— — oo 



With the choice of phase 



x may be further expressed as, 



[since Im(x n ) = 0] 



x = x n cos nq 

n— — oo 



Similarly using Eq.(^13T) and ( p4a| ) we get for v n — — V— n 

oo 

v = w n sin 



n— — oo 



Inserting Eq. fl3.11a| ) and ( |3.11b|) in Eq.flSj) and (gj) we obtain, 

OO OO »t 

J = —i nx n exp(in(f)) / drj3(t — r)iVn exp(im<? 

n— — oo m— — oo ^ 
oo 

H-zS(i) nx n exp(m</>) 



n— — oo 

oo oo 



= u)(J) + J] ^ ^jcxp(in^) / 

n— — oo m— — oo ' l; 



drf3(t — T)v m exp(imq 



dx Ti 



n— — oo 



(3.14a) 



(3.14b) 



(3.15) 



(3.16) 



where, we have expressed S(t) as a sum of two terms; the noise due to heat bath, F(t) and quantum correction term, 



S(t)=F(t)+g(t) 



(3.17) 



In the equations of motion (3.15) and (3.16), the argument of the damping mem ory k ernel (3 is (t — r ). Now (3 decays 
to zero in a time r c (the correlation time). So, to deal with the integrals of Eq.( 3.15| ) and Eq.( 3.16| ) ; it is reasonable 
to divide the range of integration into two parts: (a) \t — r\ < t c and (b) t ^> r c . Thus following CN [|10| we can write 
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<f>(t) = 0[r + (t- t)} 



(*-r) , 



neglecting higher terms of t c . It follows that, 



4>(t) ~ </>{t) - (t - r)c^ 
and v m (r) ~ u m (t) 



(3.18) 
(3.19) 



Eq.( |3.18 ) and Eq.( 3.19| ) are reasonable approximations so far as the integrals of Eq.( 3T5| ) and Eq.( 3.16| ) are concerned. 
Within the integral, we therefore manipulate the behaviour of <f> and v m for a time upto which f3(t — t) exists and 
also for the observational time at which [3 has decayed to zero. So, more specifically we can write for \t — t\ < t c , 



t rt 

dr(3(t — r)u m (r) exp[im0(r)] ~ v m (t) exp[im(f)(t)] / dr/3(t — r) exp[— im(t — t)uj] 
(i ~ Jo 



and for t ^> t c , using (|3.6|) we have 



dr/3(t - r)v m (r)exp[im^>(T)] ~ v m (t) exp[im<j)(t)}(3 m (u) 



(3.20) 



(3.21) 



Putting Eq.(3.21) which takes into account the observational time scale, in Eq.(3.15) and Eq.(3.16) we get, 

oo oo oo 

+ iS(t) nx n exp(m^) (3.22) 



j = -i 2J X! na; n w m/3m(w)exp[i(n + m) 

n— — oo m=—oo 



CO oo 



and 4> — lo(J) + x' n v m j3 m (uj) exp[i(n + rri) 



n— — oo m— — oo 



n— — oo 



where x' = 



9J 



(3.23) 
(3.24) 



Our next task is to formulate the Fokker-Planck equation. To this end we note that Lax j47j] had prescribed a method 
for deriving Markovian Fokker-Planck equation from a classical Langevin equation with short but finite correlation 
time. Although the procedure can be extended to higher order iteration scheme to include non-Markovian effects we 
adopt the method advocated by Carmeli and Nitzan |io| for their classical theory. This is based on Kramers-Moyal 
expansion of the transition probability which connects the probability distribution function P( J, <fi, t) at time t with 
that of P(J, 4>,t + t) at a later time t + r for small t, given that we know the moments of the distribution. For details 
we refer to Risken |59| The time evolution of the probability distribution P(J, <j>, t) is determined by the equation, 



dP 
~dt 



lim 



oo 
n=l 



1)" 



E 

(m,k— 0);(m+fc— n) 



d_ 

dJ 



d 



{((Aj t r(A0 t ) fc > SJ p} 



(3.25) 



where 



AJ t = AJ t (r) = J(t + t) - J{t) 
Acb t = A^(r) = <t>{t + r) - 4>(t) . 

At this juncture it is worth recalling that r is the coarse-grained time scale over which the probability distribution 
function evolves, whereas t c is the correlation time, which due to low damping is much smaller than r. The low value 
of 7 prompts us to take 7 -1 as the largest time scale for the entire problem. On the other hand, the reciprocal of 
the frequency of oscillation, i.e., w^ 1 , is the smallest time scale. Our task is therefore to evaluate the moments of the 
form ((AJ t ) m (A(f> t ) k ) s where our definition of average (. . .)s is given in (2.9). 



To evaluate the moments we make use of the following standard procedure |H 

AJ t (r) = / dsj[J{t + s),<j>(t + s),t + s] 



A0 t (r) = / ds(j)[J(t + s),(f)(t + s),(t + s) 
Jo 



(3.26) 
(3.27) 
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where the forms j and <\> are given by Eq.(3.22) and ( 3.23| ), respectively. The iterative equations are given by, 



AJ 4 (0 (r)- / dsJ[J{t) + AJ^'is), 4>(t) +Acj>l'-- x >(s),t + s 



(i-i)/ 



A#(r)= / dscj>[J(t)+AJr x '(s), <Ht) + A<l>r iJ (s),t + s 



(3.28) 
(3.29) 



where (Z) denotes the Z th iteration stage. 

The non-Markovian nature (i.e., r c is finite and r c < t ) of the present problem allows us to consider, in principle, all 
orders of r in Eq.(3.25). But, since dP/dt is evaluated in the limit r — * 0+, terms linear in r, i.e., the coarse-grained 
time scale, are taken while all the higher powers are neglected. We now introduce the following abbreviations 



c„(J) = inx n (J) 
dx n (J) 



Hn{J) 



dj 

inx n { J)v m (J)P m [uj(J)} 

dXn(J) 



dJ 



v m (J)f3 m [oj{J)} 



(3.30) 
(3.31) 
(3.32) 
(3.33) 



Substituting Eq.( p.30| )-( p.33| ) in Eq. fl3.22j ) and fl3.23| ) we obtain the quantum equations in the form of classical equations 
ofCNfHj: 



OO OC 



J = — / J B nm (J) exp[z(n + m)cj)] + S(t) ^ cr n (J) exp(m<£) 

7i — — oo m= — oo n— — oo 

oo oo oo 

and = cj(J) + C nrn (J) exp[z(n + m)<j)] — S(t) ^n(J) exp(in^) 

n— — oc 'm— — oc n— — oo 

From Eq.( |3.28 ) and (3.34) we get the explicit structure of AJ t as, 

OO OO 

AJt{r)=- £ £ / ^nm[J(t)+AJ t (s)]exp{i(n + m)[0(t) + A0 t (s)]} 

n— — oc m— — oc ^ 
oo „ T 

+ J2 dsS(s)<j n [J(t) + AJ t (s)}exp{in[<p{t) + A4> t (s)}} . 

n— — oo ^ 

Similarly from Eq.(^2S|) and ( p3E| ), A^ is given by, 

A0 t (r) = / dsuj[J{t) + AJ t (s)] 
Jo 

OO OO 

+ £ £ / rfsC nm [J(t) + AJ t (s)]exp{i (n + m)[<p(t) + A0t(s)]} 

ti=— oo m=— oo 
oo „ T 

- V / ds5(s) j u„[J(t) + AJ t (s)]exp{m[^(t) + A0 t (s)]} . 
„ Jo 



(3.34) 
(3.35) 



(3.36) 



(3.37) 



For beginning the systematic iteration procedure given by Eq.( 3.36 ) and ( 3.37 ) we initialize the zero order iteration 
stage as, 



AJ f (0) (r) =0 
and A0| O) (r) = u[J(t)]T = , 



(3.38a) 
(3.38b) 



The entire process of iteration involves cumbersome but straightforward calculations, some relevant details of which 
appear in the Appendix-A. Here we state only the main results. 

Inserting Eqs. (3.38a) and (3.38b) in the right hand side of Eqs.( 3.36| ) and (3.37) we get the results of the first order 
iteration. Thus, 
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OO OO „-7- 

AJ t (t) = — r B n _- n + cr„ exp(m</>) / dsF(s) exp(inu>s) 

Ti— — oo n— — oo ^ 

OO _ r 

+ tr„ exp(m</>) / dsg(s) exp(mws) 



(3.39) 



n— — oo 



and 



OO OO 

A^j 1 -* (r) = lot + t C n -n — J]] /i„ exp(in</>) / dsF(s) exp(mws) 

n— — oo Ti— — oo 

oo „ T 

— [i n exp(in(f)) / dsg(s) exp(inu>s) 



(3.40) 



where in writing Eq.(3.39) and (3.40) we have used Eq.(3.17). For the second iteration we put Eq.(3.39) and (3.4C) 

back into Eq. fl3.36D and ( |3.37| ) and thus obtain AJ t (2) (r) and A<^ 2) (t). Putting them back into Eq.( |3.3q ) and ( |3.3t| ) 

again we obtain AJj (r) and A0j (r) . These are p resen ted in some details in the Appendix. 

In calculating the moments as demanded by Eq.( 3.25| ), we have neglected all higher powe rs (n > 2) of r and 
l/u. The reason for doing this in case of r is clear from the limit imposed on r in Eq.( 3.25 ). For l/u> also, this 
approximation is legitimate since l/u is the shortest time scale of the problem (see inequality (|3.7|)). The final results 
for the moments are: 



([AJ t (r)] 2 ) s ^4r^n 2 K| 2 QM 



<[A<Mt)] 2 ) s 
[AJ t (r)][A^(r)]) s 
(AJ t (r)) s 



= 4 -E 



d.j 



= 



n=l 

(A^t(r))s = tJT + r n 



dJ { n dJ 



- T .fo^otc 



(3.41) 

(3.42) 
(3.43) 
(3.44) 

(3.45) 



where 



k = 


/ dt(3(t) cos(ncjt) 
Jo 


(3.46a) 


Pn = 


poo 

\ dt(3(t) sin(nwi) 


(3.46b) 




/•OO 

/ dtC(t) cos(nwi) 

./O 


(3.46c) 


<% = 


/>oo 

/ *C(t) sin(nwt) 

«/ 


(3.46d) 



and 



Also 



/o 



E 

n— 3 



5 r 



{n-l)\dt 



V n [x{t)]{{8X n -\t))-{8X n -\Q))) 



(3.46e) 



(3.46f) 
(3.46g) 



Some remarks are needed in connection with Eq.(3.41) to Eq.(3.45). Let us now examine how the quantum notion 
is implied in Eqs.(3.41-3.45). First, all the moments are the functions of the Fourier components x n (J) where J is a 



10 



quantum c-number. Second, the moments are crucially dependent on the Fourier components of quantum correlation 
function C(t) of the heat bath. In the classical limit x n (J) becomes the functions of the classical action variable J 
and also C n (ui) reduces to C n (u) — fcsT/?„(u;). We thus obtain, 



(lAJ t (T)} 2 ) s = 4Tk B TY,n 2 \x n \ 2 fa(u,) 



n=l 

oo 



([A^(r)] 2 ) s =4rfc B r^ 



dJ 



(3.47) 



in the high-temperatu re li mit Tilj <C ksT. The last term in Eq.(3.45) is due to a correction t o fre quency lo and is of 
pure quantum origin ( 3. If ). /□ is precisely the co-efficient of t in the Taylor expansion of Eq.( 3.1f ). t c is the 'cut-off' 
time upto which the quantum fluctuation remains linear in time and is approximated as ~ 1/w, as allowed by the 
time scale of the problem. Thus the quantum character of the nonlinear system enters into the description in two 
different ways. First, classical Hamiltonian gets modified by quantum corrections at t — (see Eq.( |3.3| )). This makes 
action angle variables bear quantum signature. Second, the quantum correction for t > as contained in g(t) makes 
its presence in p hase drif t term in ( 3.45 ). 

Inserting Eq.( |3.41 ) to ( 3.45| ) in Eq7(3.25) and thereby disregarding terms with n > 2 with the following definitions 



e(J)=2]>> 2 M 2 ^M 



71=1 



r(j) = 2£ 



dx r 



dJ 



O(J) = oj + ]T \ 



d 

dJ 



fis d\Xn\ 2 



(3.48) 
(3.49) 
(3.50) 



we obtain the Fokker-Planck equation for P( J, 4>,t) as, 

dP(J,<f>,t) d 



dt 



dJ 
r(j) 



2j2n 2 \x n \ 2 ^) 

n=l 

d 2 p 



dP 



+ w(J) }P 



(3.51) 



If the distribution function is initially independent of <f> it satisfies the quantum diffusion equation in action space 



dP(J, t) _ d_ 
dt ~ dJ 



(3.52) 



where by virtue of ( |3.5d| ) we write 



A = A(aJ) 



C c n (uJ) 



or A = hZJ [n(ZJ) + 1 /2] 



(3.53) 



Here lo is the linearized frequency and A plays the typical role of fegT. We have 

dH _ dE 
~dJ ~~dJ 



LO(J) 



Expressing 



we have 



lj(J) = v{E) 



(3.54) 



(3.55) 
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With this transformation the quantum Kramers' equation for energy diffusion [Eq.(3.52)] looks like, 



dP{E, t) _ d_ 
dt ~ dE 

where the diffusion coefficient is given by 



D(E) = v(E)2hU n(ZJ) + 



E" 2 ^«i 2 / 

n=l J ° 



dtj3{t) caa[nv(E)t] 



(3.56) 



(3.57) 



Eq.(3.5t) is the first key result of the present paper. The equation is valid for arbit rary temperature and noise 
correlation. The prime quantities that determine the equation for energy diffusion ( |3.56| ) are the diffusion coefficients 
D, the quantum analogue of k T, A a n d the frequ ency of the dynamical system, v(E). It is important to note that 
all the qu antit ies as defined by ( |3.57| ), ( 3.53 ) and ( 3.54 ), respectively contain quantum contributions. In the classical 
l imit Eq.(3.53) reduces to kT when n(oJ) 1/2 and n(oJ) (= [exp(ftZJ/kT) — 1] _1 ) ~ kT/huJ. Since by virtue 
( Ijg ) v{E]~= uj{J) = dH/dJ with H defined as H = {l/2)v 2 + V q {x) where V q {x) includes quantum corrections 
over the classical potential V{x) according to Eq^^), v{E) reduc es to classical frequency in the classical limit as 
ft — * 0. Although the expression for the diffusion coefficient (3.57) looks a bit complicated and formal due to the 
appearence of the Fourier coefficients x n in the summation, it is possible to read the various terms in D(E) in the 
following way. D(E) is essentially an approximate product of three terms, huJ[n(uj) + 1/2], J °° dt0(t) cos[nv(E)t] 



and K£)£~X 



where the n dependence of the latter two contributions have been separated out for 



interpretation. The integral is the Fourier transform of the memory kernel, while the sum can be shown to be equal 
to J ( Appendix D of Rcf. |l(| ), which is the quantum action variable. In the classical limit the quantum diffusion 
coefficient D(E) therefore clearly reduces to the classical diffusion coefficient of Carmeli and Nitzan p0| . A few 
further remarks on the related issues may be made at this point. To address the problem of nonequilibrium quantum 
tunneling above cross-over temperature (nonequilibrium situation arises due to the significant growth of population 
above zero levels at temperature above cross-over since the dissipation is very weak), several authors [|38|-fl3| have 
advocated the use of a probability function per unit time (of finding the system in the barrier region near a classical 
turning point with energy E) which obeys an integral equation p3fl whose the differential approximation leads to 
an equation similar (not the same) to Eq.( 3.56] ). The notable difference is in the fact that the former equation is 
applicable above cross-over temperature while Eq.(3.56) works at all temperature down to vacuum limit. The energy 
loss coefficient in ( |3.56 ), i.e., D(E)/A when put into the form [ using (3.53) ], 



D{E) 
A 



dt(3{t)v{E) 2 



n 2 \x r , 



' cos[nv{E)t] 



is comparable to that of Griff et al |43|, the integrand without (3{t) being a function of action (or equivalently energy). 
These results can be utilized as a consistency check of the present scheme. 



IV. QUANTUM ENERGY DIFFUSION CONTROLLED RATE OF ESCAPE 



The classical treatment of memory effects in the energy diffusion controlled escape is now well-documented in the 
literature p0| , p5|p^ ] . To address the correspo nding problem in the quantum domain we start by recasting the Kramers' 
equation in the energy diffusion regime |Eq.( |3.52| )1 in the form of a continuity equation to obtain, 



dP(E,t) , d jE 



= 



dt dE 

where je is the flux along the energy co-ordinate at thermal equilibrium and is given by, 



Je 



-D(E) 



d_ 1 

dE + A 



'(E)Pst{E) 



(4.1) 



(4.2) 



where P s t is the stationary probability distribution. For zero current condition, we have the equilibrium distribution, 
P eq at the source well as, 
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N- 1 

P eq {E) = —— exp(-£/A) 



HE) 
N- 1 
v{E) 



cxph(S c + AQ)/A] 



(4.3) 



where we have split the energy into classical ( E c ) and quantum ( AQ ) parts, the contribution arising from the 
latter being very small. We now define the rate of escape k as flux over population 



(4.4) 



where 

n a = total population at the source well 
P{E)dE 

Here E^ is the classical value of the activation barrier. Following BHL Es| we use a Kramers' like ansatz 



P{E)=r ] {E)P eq (E) 



to arrive at 



JE = -D(E)v{E)P eq {E) 



dr](E) 
dE 



(4.5) 



(4.6) 



(4.7) 



Integrating the above expression from E = Ei ~ A (see Eq.( |3.53" ) ) to E = E£, one derives an expression for energy 
independent current je (with E < E£) as, 



fo(A) - v(E c b )} 



dE 



A D{E)v{E)Pe q (E) 

[l- V {EZ)]D(E£ h 



A 



-S£/A 



(4.8) 



where we have used the boundary condition 7?(A) ~ 1. 

Following the original reasoning by BHL we now allow an outflow j out from each energy range E to E + dE, with 
each E satisfying the condition E > E£. Then we can write, 



d Jout = au(E)ri(E)P eq (E)dE 
which is compensated by a divergence in the vertical flow, 



dE 



-av(E)r)(E)P eq (E) 



(4.9) 



(4.10) 



Here a is a parameter which has been set approximately equal to one by BHL. Inserting the expression for non- 
equilibrium current Eq.(4.7), we obtain the ordinary differential equation for rj{E) as, 



Within small energy range above E% one can assume essentially a constant diffusion co-efficient i.e., 



dD(E) 



dE 



for E > Ei 



Substituting a trial solution of the form rj{E) = Cexp(s_E/A) for s < 0, in Eq.( 4.11 ) we have, 



4aA 2 



2 \ 1/2 



- 1 



(4.11) 



(4.12) 



(4.13) 
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Sett ing r}(E) = r)(E£) exp[s(E — Eg)/ A] and putting this into Eq.(4.7) and comparing this with the right hand side of 
Eq.(4.8) we have , 



r}(E%) = 1/(1 - s) fors<0 



Thus, escape rate k can be obtained as, 

k = je 



r 1 {E)P eq {E)dE 



Making use of (4.14) in (4.8) and the resulting expression for j in (4.15) we obtain 



k = 



1 - s 



f E » r,(E)P eq (E)dE 
(N-i/A)D(E<)eM-E c b /A) 



(4.14) 



(4.15) 



(4.16) 



For the dynamics at the bottom we have i] — > 1. Recalling that E = E c + AQ, where AQ is the quantum contribution 
to classic al en ergy, we expand £ in a Taylor series. Here E c = (v 2 /2) + V(x) and AQ(x) is the quantum correction 
terms in ( |3.3| ). Retaining terms upto the second order in x, and making harmonic approximation around the bottom 
of the well at x — 0, we get 



v 1 

E = 2 + 2 * + AQ ° 



AQ' Q x + ^AQ'qX 2 



(4.17) 



where loq corresponds to the frequency at the bottom of the classical potential V(x) at a 
d 2 V(x)/dx 2 \ x= o. The subscript zero in AQo, AQ' and AQ ' are the quantities evaluated at this point. 
Now n a , the total population at the source well, can be evaluated as, 

oo />oo 



so that 



— oo ^ — oo 



P eq (E)dxdp 

/oo />oo 
exp(-v 2 /2A)dv / exp 
-oo J — oo 



Thus 



2ttA 



AQq 



exp 



AQo 
A 



\(\ulx 2 + AQ' Q x+ l -AQ'<x 2 



d/x 



2A(lo 2 + AQ%) 



(4.18) 



So, the quantum non-Markovian rate of escape from a metastable well in the low-friction regime is given by, 



{l + (4aA 2 )AD(£g)} 1/2 -l 
{l + (4aA 2 )/D( J B 6 c )} 1 /2 + i 

AQo ' 



x exp 



DjEf) y/g + AQ>> 
A 2 

(AQo) 2 



2tt 



2{u 2 + AQq) 



(4.19) 



The above expression is the second key result of the paper. It has the form of the celebrated Arrhenius expression for 
rate coefficient with the classical activation energy E£ in the exponential factor and a complicated A and D dependent 
quantity in the pre-exponential factor. As noted earlier in the detailed discussion of quantum diffusion coefficient in 
the context of Eq.( 3.56| ), the diffusion coefficient D(E) is contributed by the three factors and has to be evaluated at 
the barrier top. The main effect of the pre-exponential factor is that the rate becomes proportional to the damping 
coefficient and the me mory kernel results in the decrease of pre- factor for increasing correlation time. The structure 
of the rate expression (4.19) suggests that it has the same form of the pre-exponential factor as that of Hanggi and 
Weiss although its content is quantum mechanical in character. The quantum mechanical content of the rate 
exp ressio n lies in se veral quantities, e.g., quantum diffusion coefficient D(E^), quantum analogue of fcgT, A as given 
by ( 3.57 ) and ( 3.53 ), respectively. The frequency at the bottom of the well, cuq as well as the classical activation energy 
Eg get modified by quantum correction AQq and AQ terms. The result of Hanggi and Weiss Q for the classical 
non-Markovian case can then be appropriately recovered. It is thus apparent that quantum correction terms in the 
exponential factor is (4.19) depends on the nature of the potential which in turn determines the rate. In what follows 
in the next section we illustrate the results with a specific cubic potential of Kramers' form. We mention in passing 
that throughout the treatment the noise intensity needs to be small for the result (4.19) to be a good description of 
the activated process controlled by energy diffusion. 
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V. AN EXAMPLE WITH CUBIC POTENTIAL 



We consider a model cubic potential of the form V(X) — -~( 1/3) AX 3 + BX 2 . A and B are two constant parameters 
of the problem with A > and B > 0. Then by virtue of Eq.(3.3) we have the c-number form of quantum potential 



V q {x) 



A 



x 3 + Bx 2 - A(SX 2 {0))x + constant 



(5.1) 



so that the time independent Hamiltonian is given by, 



H(x,v) 



A 



Bx 2 -Cx = E 



(5.2) 



where C — Ah/ (2v2~B). We have used (SX 2 (0)) = h/(2^/2B), the minimum uncertainty, and ignored the constant 
part in V q {x). Here C refers to quantum contribution to classical potential due to which the minimum, the metastable 
point corresponding to V q (x) shifts to x$ = C/(2B) (with respect to the corresponding classical metastable minimum 
at x = Q). 

Linearizing the potential V q {x) around xq we obtain 



/ AC\ 
V q (x) = V q (x ) + {B- —j (x - xo) 

We then calculate the action J, the usual form of which is given by, 

fX 2 



J = 2 



•jdx 



(5.3) 



(5.4) 



x\ and X2 are the two turning points of oscillation for which v is equal to zero and they jointly correspond to a 
particular value of the system energy E. In principle, they are the first two roots (in ascending order of magnitude) 
of the cubic equation 



.4 



Bx 2 + Cx + E = 



(5.5) 



the third root being irrelevant for the present purpose. x \ an d X2 however can be approximately calculated by simply 
putting V q (x) = E (since v — is the turning point) in ( |5.3j ) and solving the resulting equation for x. 



Xl.2 — £o T 



E-V q (x ) 
B 



1/2 



Putting the value of v from Eq.(5.2) in Eq.(5.4) we get the action integral in the form, 



J = 2V2 



(E - Bx 2 



jx 3 + Cx 



1/2 



dx 



Putting Eq.(|3.54|) in Eq.(|3.57|) we can express the quantum diffusion co-efficient in terms of the action as, 



D(J) = uj(J)2huJ 



E 



n 2 \x n \ 2 



dt(3{t) cos[nw( J)t] 



(5.6) 



(5.7) 



(5.8) 



where we have replaced D(E) by D( J) to emphasize the change made in the argument. Furthermore for unit mass 
of the Brownian particle we may write |lCfl 



;(J) ri 2 \x n (J)\ 2 = J 



n— — oc 



Putting Eq. (J5.9|) in Eq.flo.Sj) the diffusion coefficient can be approximately expressed as 

1" 



D(J) ~ 2Jhuj 



n(u> ) + 



dtp(t) cos[nu;(J)t] 



(5.9) 



(5.10) 
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For the present form of model potential we also have AQ(x) = —Cx fo r wh ich AQ g — , AQq = and AQq = — C. 
With these expressions for quantum contributions and making use of Eq.( 5.1C ) in Eq.( 4.1S ) we have the final expression 
for the escape rate as, 



2tt 



{1 + (4aA 2 )/D(J b )} 1 / 2 - 1 
{1 + (4aA 2 )/ J D( J b )}i/2 + i 



D{J b ) 



Al 4 + ^ 



(5.11) 



Here Jf, denotes the value of the action of the system at the barrier top. It should be noted that it includes both 
the classical and quantum contributions. corresponds to classical activation energy which gets modified by a 
contribution due to quantum correction entangled with the nonlinearity of the potential. It is important to note that 



the positivity of the factor (A 2 ft 2 )/ (4cjgA) in the exponential in (5.11) results in a larger effective activation barrier 
which causes a net reduction of the full rate below its corresponding classical value. This is in good agreement with 
the earlier observation by Griff et al J43| and is somewhat counterintuitive - as emphasized by Hanggi et al Q to the 
fact that full rate comprises classical rate plus zero temperature tunneling. The quantum reduction of total rate in 
the weak friction regime is a manifestation of interplay of thermal noise and quantum fluctuation and is expected 
to be pronounced for systems with flat barriers, commonly encountered in absorption-desorption processes in surface 
phenomena Q). 

VI. CONCLUSIONS 

Based on a true quantum phase space distribution function and an ensemble average procedure we have derived a 
generalized Kramers' equation for energy diffusion and analyzed the quantum transmission coefficient associated with 
the rate coefficient within a full quantum mechanical framework in the low friction regime. The present formulation 
is a complementary follow-up to our recent work [fj2| on quantum Kramers' theory in the spatial diffusion limited 
regime. The main conclusions of this study are the following: 

(i) The proposed Kramers' equation in the energy diffusion regime is an exact quantum analogue of non-Markovian 
classical Kramers' equation derived by CN |l(| in eighties. The equation retains its full validity both in the classical 
and vacuum limits at arbitrary temperature and noise correlation of the heat bath. 

(ii) The generalized quantum rate coefficient for the decay from a metastable well reduces to Kramers'-Hanggi- Weiss 
rate ]35|] in the classical limit and to pure weak dissipative tunneling rate in the quantum limit at zero temperature. 

(iii) While in the intermediate to strong damping regime the total Kramers' rate comprising classical as well as 
quantum rate is always higher than the corresponding classical rate, the notable feature in the weak friction regime ( 
for a metastable potential ) is a net quantum reduction of the total rate below its corresponding classical value. This 
is in conformity with the earlier observation in this context |S] . 

(iv) While the existing methods of calculation of quantum Kramers' rate is based on path integral techniques 
60|j6lt| , we rely on a canonical quantization procedure and true probability distribution function of c- number 
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variables. To the best of our knowledge the implementation of a differential equation and its solution as a boundary 
value problem have not been tried up to date for quantum Kramers' problem. The methodology as pursued here 
allows us to apply classical techniques for the quantum problem of barrier crossing dynamics. 

(v) The quantum effects appear in the present formulation in two different ways. The nonlinear part of the potential 
of the system gives rise to quantum dispersion, while the heat bath imparts quantum noise. An important advantage 
of the present method is that it is possible to incorporate quantum corrections to all orders and one need not invoke 
any semiclassical technique which is almost always used in the practical evaluation of the formal functional integrals. 

The present scheme of mapping of the quantum theory of Brownian motion in energy space into a classical form 
offers an opportunity to generate quantum noise as classical c-numbers and study numerically the quantum stochastic 
dynamics independent of path integral Monte Carlo techniques []62H64| . We hope to address this issue in a future 
communication. 
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APPENDIX A: CALCULATION OF MOMENTS 

Some details of the calculations involving the iterations (upto the third order) for determination of the moments 
have been shown here. The procedure followed here is that of classical theories of Carmeli and Nitzan We have 
stressed the steps for which the quantum contributions form essentially new content. 

1. First Iteration 

Inserting Eq.( |3.38gJ ) and Eq. Q3.38bp in Eq.(|3.36|) we have, 



AJj (r) = — ^2 ^2 B nm exp[i(n + m)(j)] / ds exp[i(n + m)uis] 

n— — oo m— — oo ^ 
oo ~ T oo „ r 

+ a n exp(in4>) / dsF(s) exp(inu>s) + a n exp(m^) / dsg(s) exp(mus) (Al) 

n=-oo ^° n=-oo ^° 



where we have suppressed t he ar guments of quantities B nm , a n and (f> for the sake of brevity. Now, from Eq.(3.7) and 
the argument following Eq.([3.25|) we can infer that cut 3> 1. So we can write 



ds exp[i(n + m)ujs\ ~ T<5„,_ m (A2) 

oo 

for which ^ B nm exp[i(n + m)]r(5„,_ m = r B„,-„ (A3) 



n.m— — oo n— — oo 



Thus, 



OO OO ,,7- 

AJ^^r) = — t B nt - n + cr„exp(m0) / dsF(s) exp(inujs) 

71— — OO 71 — — OO ^ 

cr„ exp(zn0) / dsg(s) exp(inuis) (A4) 
,- Jo 



n— — oo 



Similarly from Eq.(3.37) we have 

oo oo „ T 

A0( 1 ' ) (t)= w t + r ^2 C n ,- n — ^2 /i n expiiruf)) / dsF(s) exp(mws) 

n— — oo n— — oo ^ 

— /i n exp(m</>) / dsg(s) exp(mws) (A5) 

n=-oc Jo 

Along with these we will also require the statistical averages of the above expressions for calculation of the moments. 
Thus, 



OO OO „-J~ 

(AJ, (1) (t)) s = -t ^2 B n - n + ^2 cr„exp(m0) / dsg(s) exp(mws) 



(A6) 



7i— ~ oo n — — oo 



where we have used (2.1C). We can proceed further with Eq. (|A6|) and cast it in a more transparent form. In this 



context it is w orth mentioning that so far as the quantum noise g(t) is conc erned, we take only the significant terms 



from Eq.( 3.1i ) which result in terms linear in r. Referring back to Eq.( 3.32 ) we then have 



DC 



n— — oo n— — oo 



T ^2 B n,-n = -T X! n2uj \x n \ 2 (3 n (Lj) 

71— — OO 

OO 

-2rJ2 n2 "\ x n\ 2 Pn(") (A7) 



71 = 1 



17 



using Eq. fl3.46a ). The second expression of Eq.(A6) has been shown to be negligible in Appendix-B. Thus 



n=l 



(A8) 



Similarly from Eq.flAS) we have, 



(A^ 1) (r)) s = u,t + T J2nuj^^$ s n - rf^t c 



d\x n \ 2 
dJ 



(A9) 



In deriving Eq.(A9) we have used Eq.(3.31) and ( [3.33 ). Otherwise, the way lead ing to Eq.(A9) is similar to that of 
Eq. ([AsI) . The origin of the last term is the same integral that appears in Eq.(A6) and is shown in Appendix-B. 



2. Second Iteration 



Here we in sert Eq.(A4) and ( |A5| ) into the right hand side's of Eq. ( 3T36| ) and Eq.( |3.37 ). Keeping in mind Eq.(3.38a) 
and ( 3.38b ) we expand the functions of AJ t (s) and A0 t (s) , viz. B nm , C nm , cr n and keeping terms upto the first 
order only. Thus, 



B nm [J(t) + AJ t (s)} =B 

I! HI 

with B nm = B nm [J(t)] 
dB nr 



and B' = 



dJ 



(AlOa) 
(AlOb) 

(AlOc) 



j=j(t) 



The expansion has been done about AJ^°\s) = and A^'fs) = uis. The same expansion procedure is followed for 
C nm , cr„, fi n and the exponentials as well, occurring in Eq.( 3.3(f ) and (3.37). We therefore obtain, 



oo oc oo 



AJ t (2) (r) = AJ f (1) (r) " E E E \ B '™ Gl ~ l{ ~ n + m ) B nm»l] exp[z(n + m + l)<t>] 

n=— oo m— — oo l — ~ oo 

„T „ s OO OO OO 

x ds dsiF (si) exp[i(n + m)u)s + ilu>s\] — [B' nm ai — 

^ ^ TL— — OO 771— — OO £ — — OO 

i (n + m)B nm iii} exp[i(n + m + l)cj)] / ds ds%g(si) exp[i(n + rri)ujs + iluisx] 

Jo Jo 

oo co _ r „ s 

+ (a' n ai — ina n fii) exp[i(n + l)4>] / ds / dsiF(s)F(s\) exp(inuis + ilu>s\) 

r l= -oo(=-oo " ' Jo Jo 

CO CO pT pS 

(cr^c; — ina n fii) exp[i(n + l)<f>] / ds / dsiF(s)g(si)exp(inuis + ilu)Sx) 
„ Jo Jo 



n— — oo — 
oo oo 



UP OO p T pS 

(cr^cr; — ina n fii) exp[i(n + l)4>] / / dsig(s)F(si) exp(inu>s + ilu>s\) 
= -oo(=-oo ' Jo Jo 

CO CO ~ T pS 

(ct^itj — iruj n ni) exp[i(n + / ds / dsi<?(s)<?(si) exp(mws + Uujs\) 

Jo Jo 



n— — oo /= — oo 



(AH) 



The derivation this form of Eq.fAll) req uire s two i mport ant steps to be followed. The first step is just a Taylor 
expansion. While putting Eq. (|A4| ) and ( [A5| ) in Eq.( 3.36 ) and Eq.( 3.37 ), we come across several integrals which 
contain one quantity in common in the integrands. This is an expression of the form exp[ik{4>(t) + A(f>^ (s)}], where 
k is an integer. The first part, i.e. exp[ikcf)(t)], being a function of t, can be taken outside the integral while the second 
part can be dealt with as follows: 
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exp[ik{A4>P(s)}] = exp[<fc{A0 t W (a) - A0< o) ( s )}] exp^A^s)}] 

= [1 + ik{A4 1} (s) - ljs}} exp{ikujs) (A12) 

Here we have used Eg. ( 3.38b] ) and as an essential step (second) discarded the nonlinear terms. Also, since AJj '(t) 
and A<^^(r) are of order O(r), we have 

AJ t (1) (s)exp[ifcA^ 1) (s)] ~ AJ ( (1) (s)exp(ifcws) (A13) 

Guided by the same physical reasoning, we assert that all integrals of the form JJ" ds sF(s), with F(s) being finite 
as s — > 0, yield terms of the order r™ (n > 1), and are hence neglected. In a similar way, we can obtain the second 
iteration on 4>. Thus, 

oo , T „ a 



oo PT rs 

A^ 2) (t) = A^ 1} (r) + J V cr„exp(m0) / ds / dsiF(si) exp( 

n=-oo -/O J 

00 f T /«S OO OO OO 

+ w' E cr n exp(in<^) / ds I ds x g{s\) exp(mwsi) + E E E K™ "' ~ 

n— — oo n __ QO m __ 00 _ co 

z (n + m)C nm ^i^] exp[z(n + m + / ds / ds\F{s\) exp[i(n + rnjuis + zZcjsi] 

^0 ^0 

oo oo oo „ T ~ s 

+ X X [C^ m o-; - i(n + m)C nro/ n;] exp[z(n + m + Q<ft] / ds / ds 1 g(s 1 ) 

n=-oo m=-oo i = -oo 

x exp[i(n + m)ivs + ilu>s\\ 

OO OO p T p S 

— E] E] (fi^cri — inHnfii) exp[i(n + l)tf>] I ds dsxF(s)F(s\) exp(inuis + iluisi) 

7i=-ooJ=-oo JO JO 

oo oo ~ T „ s 

E] E] (/U^,* 7 ; — in^nUi) exp[z(n + l)cj>] / ds / dsiF(s)g(s±) exp(mws + iluisi) 
= -ooi=-oo ' Jo 7o 



n— — 

oo oo 



OO OO /* 7~ /* S 

— E] E] (pJ,CT| — infi n pbi) exp[i(n + l)(f>] / ds dsig(s)F(s±) exp(mws + iluisi) 

„=- oi=- 00 Jo Jo 

OO OO /.T /-S 

— E] E] (/i^cr; — in/z n /Xi) exp[i(n + / ds / dsxg(s)g(s\) exp(inujs + iluisi) 
»=-ooZ = -oo ' ' Jo Jo 

(A14) 

We next require to calculate the averages (AJ^ (r))s and (A0^ (t))^. By virtue of Eq. ( 2TTo| ) , the former yields, 

OO oo oo 

(AJ t (2) (r))s = <AJ t (1) (r)) s - E E E " ^ n + m ) B nm»i] exp[i(n + m + 00] 



n — — oc m— — oo Z= — oo 

x / eLs / dsig(si) exp[i(n + m)ujs + ildJSi 
Jo Jo 

oo oo 



"O OO /.T 

+ 2_. E/ — invn^i) exp[i(n + l)<p] / ds / dsi(F(s)F(si)) exp(inujs + iluisi) 

„=-ooZ=-oo ' *" Jo Jo 

E^ E^ (tr^ cr; — incr n fii) exp[i(n + l)<p] / ds dsig(s)g(s\) exp(inujs + ilujsi) 
„ Jo Jo 



n— — oo Z- 



(A15) 



In Eq.(Al5) we encounter three double integrals. The first and consequently the third are shown to be negligi- 
ble in Appendix-B. For the second integral, viz. J Q ds J ds\ (F(s)F(si)) s exp(mws -I- ilwsi), or J Q ds J Q dsiC(s — 
si) exp(mws -I- Uujsi), we refer to Ref. [^0|. In Appendix-C of Ref. [|l0| this type of integral of more general form have 
been evaluated. Here we state the results only. 
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rU) 



ds / dsi / ds2 ... / dsjC(s — Sj) exp[inu>s + ilusj] 
o Jo Jo Jo 



1 



0* — 1)1 v 



dwi~ 



n. — l 



(A16) 



The important physical consideration here is to neglect the terms of order C(uj)/u>. The rest of the task amounts to 
solving the above integral, also known as th e Dirichlet's condition for multiple integrals. Thus the integrals of the 
third term in the right hand side of Eq.( A15) reduce to, 



ds I dsiC(s — si) ex-p(inojs + Ulusi) = T/3i(u)5 n .-i 

!! Jo 



(A17) 



Putting Eq.(A17) back into Eq.(Al5) yields, after a little algebra, 

oo 

(AJ t (2) (r)) s = (AJ ( (1) (r)) s + 2r^; 



71=1 



2 d\x n | 

'' dj 



C c n (co) 



which, after using Eq. (|A8|) reduces to, 

oo 

(AJ t (2) (r)) s = -2r£r 

n=l 

In a similar manner we also get, 

oo 

(A^ 2) (t)) s =ct + t£' 



uPn-C S n — 



d\x n \ 
dJ 



- T fo^otc 



(A18) 



(A19) 



Next we step out for the final iteration stage, the third one. 



3. Third Iteration 



Inserting Eq.( All ) and Eq.( A14 ) in Eq.( 3.36 ) and Eq.( 3.37 ) as before we get the expansions for AJ t w (r) and 
A(j>\°' (r)which we write in the following convenient form. 

oo oo „ r 

AJ t (3) (r) = AJ t (2) (r) - £ £ B' nm exp[i(n + ro)0] / ds[AJ { t 2) (s) - AJ t (1) ( s )] 

n=— oo m=— oo 
x exp[z(n + m)6js] 

OO OO 

- X! H i(n + m)B„ m exp[i(n + m)0] / ds[A<^ 2 } (s) - A^ (1) (s)] 

n— — oom- — oo ^ 

x exp[i(n + m)u;s] 

oo „ T 

+ <exp(m0) / ds[AJ { t 2 \s) ~ A4 1] (s)}F(s)exp(inujs) 

oo -, r 

+ ina n ex.p(in(j)) / <is[A(/)| 2 ' ) (s) — A0^ (s)]i ? (s) exp[i?icijs] 

n=-oo 

oo „ T 

+ cr^exp(in0) / ds[AJj (s) — Ajj: 1 \s)]g(s) exp(inu>s) 

n=~oo Jo 

°° PT 

mcr„exp(m0) / ds[A0( 2 ^(s) — A(/)j 1 ^(s)](7(s) exp[wws] 
Jo 



(A20) 



n— — oo 



and 
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A4 3) (r) = a4 2 \t) +J f ds[AJ t (2) ( s ) - AJ ( (1) ( S )] 

Jo 

OO OO 

Yl C ex P[ ! (« + m )^l / ds[AJ t (2) (s)- AJ ( (1) (s)]exp[i(n + m)ws] 

t — rv» 171 — rv% 



n— — OO 771 — — OO 
OO OO 



+ Y Yl *(" + m )Cmnexp[i(n + m)0] / ds[A^ 2 \s) - A<^ x) (s)] 

n— — oom — — cjo ^ 

x exp[i(n + m)ws] 

— ^ fi' n exp(in<f>) / ds[AJ ( (2) (s) - AJ t (1) (s)]F(s) exp(mws) 

n=-oo ^° 

— Y^ in(j, n exp(in<fi) / ds[A<f)\ 2 \s) — A<j>t (s)]F(s) exp(inujs) 
n=-oo ^° 

OO - T 

— Y^ rfn exp(m^) / eLs[AJ t ^(s) — A ^ (s)] <?(s) exp(inws) 

n=-oo ^° 

m/i„ exp(m(/)) / <is[A(/)| 2 ' ) (s) — A</>j (s)]g(s) exp(inws) 
—~ Jo 



(A21) 



n— — oo 



It is needless to carry out further iterations because they yield terms of order C(uj)/uj even in the order of r, and hence 
are negligible. From this point we therefore proceed to calculate the moments ((A J) m (A(f>) k ) s as the right hand side 



of Eq.(3.25) demands. While performing the averaging, some integrals occur in common with both (AJ t (r))s and 
(A(f) t (t))s- They are as follows:- 



h 
h 
h 
h 
h 
h 
It 



ds( 
ds( 
ds( 
ds( 
ds( 
ds( 
ds( 



Ajf( S 



Ajf( S 



A^ 2) ( s 



Ajf( s 
A^ 2) ( s 
Ajf( S 
A^ 2) (s 



AJ t (1) (.s)]>5 

AJf {s)])s exp[i(n + m)ujs] 
Acjyp {s)])s exp[z(n + m)tus] 
AJ 4 (1) (s)] J F 1 (s)) s exp(mc i ;s) 
A^ 1) (s)]F(s)) s exp(mws) 
A 1 - > (s)])sg(s) exp(mws) 
A 0t (s) exp(inws) 



(A22a) 
(A22b) 
(A22c) 
(A22d) 
(A22e) 
(A22f) 
(A22g) 



Integrals I2 to I7 are in common with (AJ t lJ '(r)) s and (A$ d; (r)) s while h occurs in the latter only. The integrands 
containing F{.) alone or in a product with g(.), on averaging gives zero by virtue of Eq.( 2.10| ). And the integrands 
containing g{.) alone or two <?'s in product, amounts to negligible contributions (see Appendix-B). T he in tegrands 
containing the correlation average of two F's and not obeying (or can not be recast) the form of Eq.(A16) have all 
been shown to be negligible by CN in their Appendix-C p0|| . Such integrals as stated above constitute the bodies 
of Ii to 7 3 . For I4 through I-j, some integrals contain product of three terms viz. (F(.)F(. )F(.) )s, (F(.)) sg{-)g(-), 
(F(.)F(.))sg(-), g(.)g(.)g(.) etc. The first two type of integrals vanish, the former due to Eq.( 2.10|) and the latter due 
to the Gaussian property of F. The last two types vanish, the latter following from Appendix-B and the former has 
been shown to be negligible in Appendix-C. The only non-zero contribution comes from ^5, the integral as a whole 

appears as uj' 07 exp(iZ0) JJ" ds J Q S dsi J Q S1 ds2C{s — S2) exp(inuis + UU0S2) , where we have used Eq.(3.5b). Using 

the result of Eq. (|Al6| ) followed by a little bit of algebra we eventually arrive at the value of I 5 . Thus, 



ITLO 



dC- n (u) 



and 



h = ~- — o-_ n exp(-m#, 

n duj 

h = h = h = h = h = h = 



(A23) 
(A24) 
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Thus we have, 



(Ajf , (r)) s = (AJ t (2) (r)) s + ]T ina n exp(in<l>)[h] 



(A25) 



n— — oo 



which, by changing the dummy index from —n to n for this symmetric summation (—00 to + 00), and employing 
Eq. (|3.12D , ( |3.30| ) and ( |3.46g|) and finally putting the value of (AJ t (2) (t)) s from Eq.( |Al8[ ) we reach at, 



(A4 3 \r)) s = -2rJ2n 2 

n=l 

Similar calculations establish the result, 



u>\x n \ 2 PU")-Jj{K\ 2 C c n (u;)} 



(A</»l 3) (r)) s = W r + r^n 



n=l 



z s d\x n \< 



dJ dJ \ 71 dJ 



C 



(A26) 



(A27) 



For the calculations of the second moments ((A Jt) 2 )s, ((A</>*) 2 )s and (A J t A(j> t )s we refer the reader to the Appendices 
B and C of Ref. p0| , where it has been clearly shown that terms which do not appear in AJ^ (r) and A^ 1 -* (r) lead 
to cross terms or square terms containing three or more integrals over one or two C functions. All of them contain 
higher powers of r or terms of order [C(uj)/lo\ n with N > 1, and are hence discarded. Thus we obtain the second 
moments as, 



((AJ t ) 2 >5 



71=1 

(AJ t A<p t ) s = 



4 -E 



dXr, 



dJ 



(A28) 
(A29) 
(A30) 



In obtaining the results of Eq.( |A28| ) to ( A30| ) we have used Eq.( AlOa) to Eq.( AlOc) and the like, along with the fact 
that higher order terms r™ and [C{oj)/u)] n (with n > 1) are negligibly small. 



APPENDIX B: EVALUATION OF INTEGRALS INVOLVING QUANTUM CORRECTIONS 

Here we calculate the terms involving single or multiple integrals of the quantum fluctuation g(.) keeping in mind 
the time scale of energy diffusion. First of all we recall the structure of g(t) from Eq.(3.U) and then refer to Eq. (|A(|) . 
On the right hand side of that equation we have an expression of the form, X^r^L-oo Un exp(m</>) /J" dsg(s) exp(mws). 
For brevity, we set 



1 



(n - l)!m! dt m 



V 



(x(t)) {(SX^it)) - (^("-^(O)))! = T mn (0) 
l. j J t=o 



so that, Eq.(3.1i) appears as 



00 00 



■9W = -EE <mT "-( ) 



m—0 n—3 



and we can write 



<j n exp(zn0) / dsg(s) exp(mcjs) = er„ cxp(zn0) 

i=-oo n= — 00 1=3 

00 „ T 
x ^^Tfcj(O) / ds s k exp(ijiws) 



(Bl) 



(B2) 



(B3) 



From Eq.(Bl) it is obvious that 
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Ton(O) = 



(B4) 



and thus it follows that the right hand side 

of Eq.( p33| ) assumes the form Xl^L-oo Sz^3 a n exp(in</>) YlkLi ^fci(O) /J" <iss fe exp(mu;s) where the summation over 
the index fe now extends from one to infinity. Within the integral the exponential part fluctuates rapidly ( lo being 
very large ), remaining finite in the limit s — + 0. As it stands even for k = 2 such integrals yield terms of order t p ( 
where p > 1 ) Jl(|. Thus, integrals of the above form are negligibly small and are hence discarded. The case of k = 1 
calls for a special attention. The above expression takes a more simple form. Thus, 

oc oo oo -, r oo oc 

cr n exp(w0) Tu (0) / dss k exp(inws) ~ ^ ^ g n exp(w^)Ti; (0) 

n=-oo 1=3 k=l n=~oo 1=3 



n= — oc Z=3 

x / dss exp(mws) 



(B5) 
Jo 

In comparison to the rapidly varying exponential part, the linear part arising from the qua ntum fluctuation does not 
alter significantly within the range of integration and hence the right hand side of Eq.(B5) can be written as 



oo oo 



[right hand side of Eq.(B5)] = a n exp(in^i)Tu(0)t c T6 n fi 

n=—oo 1=3 
oo oc 

= ^2( inx n)exp(in(/))Tii(0)t c T5 rh o [using 



n=— oo /— 3 

= o 



Here i c is the cut-off time (~ as h as been mentioned earlier in the discussions preceding Eq.(3. 

after averaging the last term of Eq.(A5) becomes, 



(B6) 
Similarly, 



oo oc 



[last term of Eq.(A5)] = ^ exp(m0)Ti;(O)t c rJ n] o 

n— — oo Z=3 
= -HotcTfo 



(B7) 



where we have used Eq.flOll) and Eq. flOfol ). 

Now consider the last integral of the right hand side of Eq.( A15 ). It is of the form JJ" ds dsig(s)g(si) exp(mws + 
iluisi). Using Eq.(B2), it can be cast in the following form:- 



ds / dsig(s)g(si) exjp(imus + Ulvsi) 
Jo 

dsg(s) exp(inujs) / dsig(si) exp(iluisi) 



oo oo 



exp(inws) / ds 



oo oo 



- E*rTmxfa(o) 

mi— O/ci— 3 



m=0 fe=3 

OO OO OO OC „ T p S 

EEE E Tlt (°) Tll '( )/ dss m exp(inujs) ds^ 1 exp(iZwsi) 

OO OO „ s 

y^ T 1 fe(0)T 1 fc 1 (0) / dssexp(inuis) / ds^i exp(^ws 1 ) [using Eq.(B5)] 
k=3k 1= 3 J ° Jo 

± ± T lfc (0)T lfcl (0) r-^(««-) _ 



exp(i/u;5i) 



fc=3 fei —3 



(B8) 



where evaluating the third bracket of the last expression we have neglected terms of order 0(l/u> 2 ). The last expression 
yields terms of negligible contribution due to the reasons given at the onset of Eq.(B5). Integrals involving three g 
functions can be similarly shown to be negligibly small. 
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APPENDIX C: TREATMENT OF INTEGRANDS OF h TO h OF APPENDIX-A 



In the discussions that appeared in between Eq.(A22g) and Eq.(A23) we mentioned of four ty pes of integrands 
which appear in the calculations involving the evaluations of the integrals I4 to I7 ( Eqs.(A22d)- (A22g) ). Let us 
recall their forms again. They are, 

(1) {F(.)F(.)F(.)) S , (2){F(.)) s g(.)g(-l (3) {F(.)F(.)) s g(.), (4) g(.)g(-U) 

It was reasoned there that types (1), (2) and (4) yield terms of negligible contributions (see Appendix-B). Here we 
clarify the way of dealing with type (3). The specific forms of integrals involving this type of integrands are:- 

ds 2 (F(s)F(si)}sg(s 2 ) exp[iw(js + nsi + ls 2 )] 



ds 2 (F(s)F(s 2 )}sg(s 1 ) exp[iu(js + nsi + ls 2 )} 



(A) 
(B) 

These two integrals occur in I4. 



ds / dsi 
Jo 

r 

ds / ds\ 



(C) 



ds / dsi 



ds 2 (F(s)F(s 2 )) s g(s 2 ) exp[iu(js + ns 1 + ls 2 )} 



This, with type (A) above, occur in 1$. And lastly the integral, 



dsx / ds 2 (F(si)F(s 2 )) s g(s)exp[iu}(js + nsi + ls 2 )] 



(D) 

occurs in Iq and I7. 

All the integrals from (^4) to (D) lead to negligible contributions. We establish this b y sho wing the case of, say, 
type (B). For doing this we invoke the quantum fluctuation-dissipation relation from Eq.( [3.5b ) as, 



(F(s)F(s 2 ))s 



1 

2 Jo 



with 



G(w') = -K(oj')p(u)')hui' 



dio' K^^p^uj^huj' coth ■ 
du)'G(tu') coso/(s — s 2 ) 

tiuj' 



Two' 
2k^f 



cos lo'(s — s 2 ) 



coth 



2k B T 



where, a superscript ""' has been added as a superscript to co to denote the frequency of the bath modes. 

1 r°° 

(F(s)F(s 2 )) s = - / eL/G(u/)[exp{iu/(s - s 2 )} + exp{-iu;'{s - s 2 )}} 
z Jo 



When Eq.(|Cl|) is put in type (B) integral along with Eq.(B2), we get the value of the integral as, 

2lu' 



p OO 

[type B]=fU du>'G{u') 
Jo 



(CI) 



(C2) 



This being of higher order in 1/u, can be discarded. It is easy to show that the other three types of integrals also 
lead to negligible contributions. The procedure is the same as that adopted above in the case of type (B). 
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